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1. Introduction 

Recently, there is a growing interest in the study of 2D gapped Dirac materials 
isoestructural with graphene. One of these materials is silicene, which is a two 
dimensional crystal of silicon, with a relevant intrinsic spin-orbit coupling (as compared 
to graphene), studied theoretically [HE] and experimentally (3]IHElE]E] • Other gapped 
Dirac materials are germanene, stannene and Pb [8|. For these systems, the low energy 
electronic properties can be described by a Dirac Hamiltonian, like in graphene, but 
the electrons are massive due to the relative large spin-orbit coupling Ago- In fact, 
we will consider the application of a perpendicular electric held Sz = Az/l {I is the 
inter-lattice distance of the buckled honeycomb structure) to the material sheet, which 
generates a tunable band gap |Ag^| = |(A^ — s^Aso)/2| (s and ^ denote spin and valley, 
respectively). There is a topological phase transition [9] from a topological insulator (TI, 
lAj-l < Ago) to a band insulator (BI, |A^| > Ago), at a charge neutrality point (CNP) 
aI^^ = s^Ago, where there is a gap cancellation between the perpendicular electric held 
and the spin-orbit coupling, thus exhibiting a semimetal behavior. 

A 2D topological insulator was theoretically studied in [10] and hrst discovered 
experimentally in HgTe quantum wells in [n]. A TI-BI transition is characterized by a 
band inversion with a level crossing at some critical value of a control parameter (electric 
held, quantum well thickness, etc). Recently, we have found that electron-hole Wehrl 
entropies of a quantum state in a coherent-state representation provide a useful tool to 
identify TI-BI phase transitions na. 

In this work we explore the connection between the TI-BI transitions in some 2D 
Dirac materials with the Shannon information entropies of the wave packet probability 
densities in position and momentum spaces and, in particular, with the entropic 
uncertainty relation. The uncertainty principle can be quantihed in terms of the 
usual variance-based uncertainty relation, AxAp > |, or alternatively by means of the 
entropic uncertainty relation [laiiiiis] which has been shown to be more appropriate 
in diherent physical situations [161 EH EH EHl EHl EH EH EH ES] . If we dehne the position 
and momentum densities of a state T as p(r) = |d'(r)p and 7 (p) = |<I)(p)p, respectively, 
with \l/(r) the position and $(p) the momentum wave packets, the entropic uncertainty 
relation is given by 

Rp-|-> D(1-|-In tt) (1) 

where D is the dimension of the position and momentum space and where Sf = 
— f /(x) ln(/(x))dx is the so called Shannon information entropy of a density /. 
The equality is reached when the wave packets in position and momentum spaces 
are Gaussians. The Shannon information entropy measures the uncertainty in the 
localization of the wave packet in position or momentum spaces, so that the higher 
the Shannon entropy is, the smaller the localization of the wave packet is; and the 
smaller the entropy is, the more concentrated the wave function is. Besides we will 
consider the relative or Kullback-Leibler entropy to characterize the topological-band 
insulator transitions (see section 0]) in these materials. 
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Ago (meV) 

1{A) 

Vp (lO^m/s) 

Si 

4.2 

0.22 

4.2 

Ge 

11.8 

0.34 

8.8 

Sn 

36.0 

0.42 

9.7 

Pb 

207.3 

0.44 

~ 


Table 1. Approximate values of model parameters Ago (spin-orbit coupling), I 
(interlattice distance) and vp (Fermi velocity) for two dimensional Si, Ge, Sn and 
Pb sheets. These data have been obtain from first-principles computations in [5] (Ago 
and /) and [Ml [50] Uf)- 


The paper is organized as follows. Firstly, in Section [2l we shall introduce the 
low energy Hamiltonian for some 2D Dirac materials (namely: silicene, germanene, 
stantene,...). Then, in Section [31 we will characterize topological-band insulator 
transitions in silicene in terms of the entropic uncertainty relation. In Section 01 the 
connection between the relative entropy and the topological-band insulator transitions 
is studied. In Section [5] we use the Heisenberg uncertainty relation to characterize 
the TI-BI phase transition. Finally, some concluding remarks will be given in the last 
Section. 


2. Low energy Hamiltonian 


Let us consider a monolayer silicene him with external magnetic B and electric 8^ helds 
applied perpendicular to the silicene plane. The low energy effective Hamiltonian in the 
vicinity of the Dirac point is given by [9] 

Hi = vf{(t^Px - ^(TyPy) - T (2) 

where ^ corresponds to the inequivalent corners K {^ = 1) and K' (^ = —1) of the 
Brillouin zone, respectively, Uj are the usual Pauli matrices, vp is the Fermi velocity of 
the Dirac fermions (see Tabled] for theoretical estimations in Si, Ge and Sn), spin up and 
down values are represented by s = ±1, respectively, and Ago is the band gap induced 
by intrinsic spin-orbit interaction, which provides a mass to the Dirac fermions. We 
are considering the application of a constant electric held 8^ which creates a potential 
diherence A^ = 18^ between sub-lattices. The value / appears in table [T] for diherent 
materials. The values of the spin-orbit energy gap induced by the intrinsic spin-orbit 
coupling has been theoretically estimated |26l [23 [28l [8] for diherent 2D Dirac materials 
that we show in Table dl 

The eigenvalue problem can be easily solved. Using the Landau gauge, A = 
{0,Bx,0), the corresponding eigenvalues and eigenvectors for the K and K' points are 
given by P EH d2] 





n ^ 0, 


( 3 ) 
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n)si 


-^+) A 


where we denote by = (1 ± ^)/2, the Landau level index n 
cyclotron frequency u = vp \/2eB/h, the lowest band gap = 
the constants and are given by EB 


(4) 

= 0, ±1, ±2,..the 
(A^ - s^Aso)/2 and 


-4? = 1 

1 ^gn(n)\/ 

[ n = 0, 

1 

f / |£^n^|-sgn(n)Asj 

B,f = 

V 2iKh 

1 

[ ^+, n = 0, 


n ^ 0, 


, n 7^ 0, 


(5) 


The vector ||?7,|) denotes an orthonormal Fock state of the harmonic oscillator. 

We will do all the numerical analysis in silicene but the results will be valid just 
replacing the corresponding parameters in Table [T] 

As already stated, there is a prediction (see e.g. [26l [23 [28l [32] ) that when the gap 
[A^^l vanishes at the CNF Ai°\ silicene undergoes a phase transition from a topological 
insulator (TI, |A^| < Ago) to a band insulator (BI, |A^| > Ago). This topological phase 
transition entails an energy band inversion. Indeed, in Figure [T] we show the low energy 
spectra [3] as a function of the external electric potential A^ for B = 0.05 T. One can 
see that there is a band inversion for the n = 0 Landau level (either for spin up and 
down) at both valleys. The energies and Eq^’^ have the same sign in the BI phase 
and different sign in the TI phase, thus distinguishing both regimes. We will observe 
a similar “inversion” behavior in entropic and variance-based uncertainty relations for 
the Hamiltonian eigenstates [dj thus providing a quantum-information characterization 
of the topological phase transition. 


3. Entropic uncertainty relation and topological phase transition 

In order to compute Shannon entropies, firstly we have to write the Hamiltonian 
eigenstates [H in position and momentum representation. We know that Fock (number) 
states \n) can be written in position and momentum representation as 

{x\n) = e {y/ux) (6) 

^2^n\y/7r 

{p\n) = -^====e-P^^‘^‘^Hn (p/\/w) (7) 

where Hn{x) are the Hermite polynomials of degree n. We will introduce the number- 
state densities in position and momentum spaces as Pn{x) = |(x|n)p and 'yn{x) = 
|(p|n)p, which are normalized according to J pn{x)dx = 1 and J ■jn{x)dx = 1. Now, 
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Figure 1. Low energy spectra of silicene as a function of the external electric potential 
for B = 0.05 T. Landau levels n = ±1, ±2 and ±3, at valley = 1, are represented 
by dashed (electrons) and solid (holes) thin lines, black for s = — 1 and red for s = 1 
(for the other valley we simply have The lowest Landau level n = 0 

is represented by thick lines at both valleys: solid at ^ = 1 and dashed at ^ = — 1. 
Vertical blue dotted grid lines indicate the CNPs separating BI (|A 2 | > Ago) from TI 
( lA^I < Ago) phases. 


taking into account Eq. (|1]), the position and momentum densities for the Hamiltonian 


eigenvectors |4] are given, respectively, by 

= ««)^Kx||n| - + (a?)"Kx||n| - (8) 

XHp) = (^?)TpII"I - + (a?)Tpll"l - P) 

We will study the position and momentum entropies 

/ OO 

p^J(x)ln(p^J(x))dx (10) 

-OO 

/ OO 

7nHp)in(7nHp))dP- ( 11 ) 

-OO 


If we make a change of variable, it is straightforward to see that S = S — In (cj). 

In £gure[2]we plot S' *5 + S' s? as a function of the external electric potential for 
the Landau levels: n = ±1, ±2 and ±3 (electrons in blue and holes in red), with spin 
up (dotted lines) and down (solid lines) and magnetic field B = O.OIT at valley = 1- 
Electron-hole entropy curves cross at the CNP lA^I = Ago- For electrons (resp. holes), 
the asymptotic entropies in position space are given by S'p|^|_j (resp. S'p|^|) for A 2 —)■ oo 
and S'pi^i (resp. S'p|^|_J for A^ —)■ — 00 . In momentum space the behavior is analogous. 
The uncertainty relation has the same value for spin up (resp. down) electrons and 
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Figure 2. Entropic uncertainty relation Sp^i + as a function of the electric 
potential for the Landau levels: n = ±1, ±2 and ±3 (electrons in blue and holes in 
red), with spin up (dotted lines) and down (solid lines) and magnetic field B = O.OIT 
at valley ^ = 1. Electron and hole entropy curves cross at the critical value of the 
electric potential = —sAgo (vertical black dotted grid hires indicate this CNPs), 
^P\r,\ + ‘5'7|n| are the asymptotic values for A^ —±oo (horizontal black dotted grid 
lines). 


A A 

^pSl SyS^l 



Figure 3. Combined Shannon entropies IT^ and IT^ of electron plus holes for n = 1, in 
position and momentum spaces, S'psi (left) and (right), respectively, as a function 
of the electric potential A^, for B = O.OIT, with spin up (blue) and down (red) and 
valley ^ = 1. The combined entropies have a maximum at the critical value of the 
electric potential Ai°^ = — sAgo (vertical dashed grid lines indicate these CNPs). 
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/KLp^i+/KLy^i 



Figure 4. Kullback-Leibler uncertainty relation IxLp^^ + as a function of the 

electric potential for the Landau levels: n = ±1,±2 and ±3 (electrons in blue and 
holes in red), with spin up (dotted lines) and down (solid lines) and magnetic field 
B = O.OIT at valley ^ = 1. Electron and hole entropy curves cross at the critical value 
of the electric potential = —sAgo (vertical black dotted grid lines indicate this 
CNPs), 


holes at the CNP point = Ago (resp. A^ = —Ago). Moreover, for each n, note that 
in the BI phase the electrons (resp. holes) uncertainty goes to greater value than holes 
(resp. electrons) uncertainty for A^ < —Ago (resp. A^ > Ago). We have checked that 
the smaller the magnetic held strength, the sharper this effect is. 

In hgure |3] we have plotted the combined entropy of electrons plus holes 


SgSl - SpBl + S, 


P-r, 


In In '—n 


( 12 ) 

(13) 


in position and momentum representation for n = 1. We can observe that the combined 
entropies exhibit a maximum at the CNPs in both representation spaces. This is a 
common feature for general n. 


4. Kullback-Leibler entropy 

The relative or Kullback-Leibler entropy is a measure for the deviation of a density /(r) 
from a reference density ^'(r) [33] is dehned as 

IklU, 9) = j /(r) In (14) 

Recently, the relative Renyi and Kullback-Leibler entropies have been found to be an 
excellent marker of a quantum phase transition in the Dicke [Ml 135] and in the U (3) 












Uncertainty relations and topological-band insulator transitions 


vibron model [SB]. In this section we will explore the utility of the relative entropy as 
an indicator of a topological phase transition. For this purpose, we will consider as 
reference densities the position and momentum densities po{x) and 7o(p), respectively, 
which are the densities for minimum uncertainty in relation ([T]). Therefore, we will 
analyze how different is a density from the minimum uncertainty density. 

We will study the position and momentum entropies 

Again, it is straightforward that ~ 

In hgure |4] we plot the sum Ixipp + as a function of the external electric 

potential for the Landau levels n = ±1,±2 and ±3 (electrons in blue and holes 
in red), with spin up (dots lines) and down (solid lines). The hgure corresponds to 
a magnetic held B = O.OIT. Electron-hole relative entropy curves cross at the CNF 
|A^| = Ago at which they reach the values ~ 0.37, 1.96 and 3.69 for \n\ = 1, 2 
and 3 respectively. For electrons (resp. holes) the asymptotic entropies in position 
space are given by Iklp\„\_i (resp. Iklp\^\) for A^ cx) and Iklp\^\ (resp. Iklp\„\_i) 
for A^ —)■ — cx). In momentum space the behavior is analogous. Note the electron-hole 
entropy inversion phenomenon anticipated at the end of Section [2l Indeed, the quantity 
^KLpp + — In has the same sign for spin up and down electrons (idem for holes) 

in the BI phase (|A^| > Ago) and diherent sign in the TI phase (|A^| < Ago). 

5. Heisenberg uncertainty relation and topological phase transition 

Entropic uncertainty relation provides a rehned version of the Heisenberg uncertainty 
relation (see [3^ and references therein); 

1 1 

AxAp > - exp [Fp -|- S'.y — 1 — In vr] > -. (17) 

It gives a stronger bound for the variance product than the standard In this section, 
we shall explore the more usual Heisenberg uncertainty relation. 

Introducing position and momentum operators through a bosonic mode, [a, a^] = 1, 
as usual: 

^ + a), P = *^ 1 ( 0 ^-a), (18) 

we can easily compute the expectation values of X and P and their fluctuations in an 
energy eigenstate 0] as 

{n\X\n)s^ = 0 , {n\P\n)s^ = 0 , 

+ 1 ) = 

Lu Z Lu 


(19) 

( 20 ) 
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Figure 5. Expectation value of the number operator N = a in the energy eigenstate 
\n)s^ as a function of the electric potential for the Landau levels: n = ±1, ±2 and 
±3 (electrons in blue and holes in red), with spin up (dotted lines) and down (solid 
lines) and magnetic field B = O.OIT at valley f = 1. Mean number curves cross at the 
critical value of the electric potential = — sAgo (vertical black dotted grid lines 
indicate these CNPs), for which = |n|/2 (horizontal black dotted grid lines) 


where 

/V*« = K«)^(|n|-l) + (B?)^|n| (21) 

is the expectation value of the number operator N = a^a in the energy eigenstate \n)s^. 
Therefore, the product of standard deviations, + is written in terms 

of solely. In Figure |5] we represent as a function of the electric potential for 
several Landau levels. We observe the same electron-hole curve inversion phenomenon 
as for the entropy curves. The quantity has the same sign for spin up and 

down electrons (idem for holes) in the BI phase (|A^| > Ago) and different sign in the 
TI phase (|A^| < A^o). 

6. Conclusions 

We have explored how different quantihcations of the uncertainty relation characterize 
a topological phase transition in a group of 2D Dirac gapped materials (monolayer 
sheets of Si, Ge, Sn, and Pb). Firstly we have inspected the entropic uncertainty 
relation. We have found that the electron-hole entropic curves cross at the charge 
neutrality poin (CNF), that is, the electron and holes with spin up (down) have the 
same uncertainty at the critical point A^ = Ago (A^ = —Ago). The combined entropy of 
electrons plus holes shows a maximum at the critical points in position and momentum 
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spaces; therefore, the combined electron plus hole density, for spin up (resp. down), 
is more delocalized at the CNP = Ago (resp. A^ = —Ago). Furthermore, we have 
analyzed the relative (or Kullback-Leibler) entropy, hnding that the sum of the relative 
entropies in position and momentum spaces identihes each phase (band and topological 
insulator) at the CNPs. For completeness we have considered the product of standard 
deviations observing the same electron-hole uncertainty curve inversion phenomenon. 
Summarizing, we have related the uncertainty principle and the topological phase 
transitions in this model. We expect that this analysis might be applicable to other 
problems to obtain a general/deepest connection between both concepts: uncertainty 
principle and topological phase transitions. 
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